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The heat capacity of a 2H-NbS2 single crystal has been measured by a highly sensitive ac technique 
down to 0.6 K and in magnetic fields up to 14 T. At very low temperatures data show excitations 
over an energy gap (2As/fcsT c ~ 2.1) much smaller than the BCS value. The overall temperature 
dependence of the electronic specific heat C e can be explained either by the existence of a strongly 
anisotropic single-energy gap or within a two-gap scenario with the large gap about twice bigger 
than the small one. The field dependence of the Sommerfeld coefficient 7 shows a strong curvature 
for both principal-field orientations, parallel (i?||c) and perpendicular (H _L c) to the c axis of the 
crystal, resulting in a magnetic field dependence of the superconducting anisotropy. These features 
are discussed in comparison to the case of MgE>2 and to the data obtained by scanning-tunneling 
spectroscopy. We conclude that the two-gap scenario better describes the gap structure of NbS2 
than the anisotropic s-wave model. 



I. INTRODUCTION 

The old concept of multiband/multigap superconduc- 
tivity [1| has found its strong experimental evidence only 
recently in the rich physics of magnesium diboride [2j. 
Consequently more superconductors are (re)considered 
along this line. One of the important examples are the 
iron pnictides, a new class of high-T c superconductors 
3 for which multigap superconductivity is suggested to 
lead to an exotic pairing mechanism with a sign reversal 
of the order parameter between separated Fermi-surface 
sheets. A revision in dichalcogenides brings more and 
more signatures of a distribution of superconducting en- 
ergy gaps, which can be either due to different gaps on 
different Fermi-surface sheets or anisotropic single gap. 

Transition-metal dichalcogenides 2H-MX2 (M = Nb, 
Ta, X = S, Se) are materials with layered structure. 
Nb or Ta atoms are trigonally prismatic coordinated by 
chalcogen atoms and metallic layers arc held together by 
weak van der Waals forces. Because of this layered struc- 
ture, electrical, magnetic, and optical properties show a 
high degree of anisotropy. NbS2 is the only member of 
the 2H-MX2 family, which does not undergo a charge- 
density wave transition [4] . This could be a reason for its 
strong anisotropy, much larger than that of NbSe2 . 

NbSe2 had been considered for a long time as being a 
conventional type II superconductor [j| . Later on, effects 
of the anisotropic and strong coupling interactions were 
taken into account [6, 7]. Recent measurements sensitive 
to the order parameter show evidence that more than 
one energy scale is necessary to account for establish- 
ing superconductivity NbS2 was also originally 



considered as just another anisotropic superconductor 
and its unusual specific heat dependence was not inter- 
preted in detail [13l . flij ]. An important breakthrough 
came with scanning tunneling microscopy /spectroscopy 
(STM) measurements by Guillamon et al. [Hj], show- 
ing strong indications for two superconducting energy 
gaps instead of a single anisotropic one. Since STM is 
a surface probe, this strong statement certainly needs 
independent support showing that the two gaps are re- 
flecting the bulk properties of the system. In this paper 
we address this issue with bulk thermodynamic measure- 
ments of the specific heat at temperatures down to 0.6 K 
and in magnetic fields up to 14 Tesla via ac-calorimetry 
technique. We find that the electronic specific heat C e 
cannot be described by the standard BCS model with a 
single isotropic energy gap. First, at the lowest tempera- 
tures the data shows that quasiparticles are excited over 
an energy gap much smaller than the BCS weak cou- 
pling limit. The overall temperature dependence of C e 
can be described only if two gaps or an anisotropic one 
gap case is considered. Second, the field dependence of 
the Sommerfeld coefficient 7 shows a strong curvature in 
striking similarity with that of NbSe2 and MgB2. How- 
ever, the anisotropy of 7 decreases with magnetic field in 
an opposite manner compared to the latter system. Fi- 
nally, the two gap scenario is supported by the absence 
of in-plane gap anisotropy in recent STM imaging of the 
vortex lattice in NbS2 [15|, and by the fact that the nu- 
merical values of the two gaps obtained from fitting our 
data, 2A s /k B T c » 2.1 and 2A L /k B T c « 4.6, are also in 
a very good agreement with the STM data. 
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II. EXPERIMENT 

Details of the synthesis of the single crystalline samples 
can be found elsewhere [lfjj]. The crystals used for the 
specific heat measurements come from the same batch 
as those used in the previous STM studies [l5|. Their 
chemical composition was checked using an energy dis- 
persive X-ray spectroscopy (EDS). Single crystals were 
also confirmed to be of 2H polytype by X-ray diffraction 
measurements. In our experiment a thin crystal with a 
well defined hexagonal shape and dimensions 500x500x30 
/im was chosen. 

Specific heat measurements have been performed using 
an ac technique as described elsewhere [l7j . The high 
sensitivity of this technique is not only very well adapted 
to measure the specific heat of very small samples but 
also to carry continuous measurements during temper- 
ature or magnetic field sweeps. We were thus able to 
obtain the field dependence of the electronic part oiC/T 
at T Pi 0.6 K which only differs from its zero tempera- 
ture limit, the Sommerfeld coefficient 7 , by about 2%. 
Measurements were performed with the magnetic field 
aligned along the two main crystallographic orientations, 
i.e. parallel and perpendicular to the basal ab plane of 
the sample. The temperature oscillations of the sample 
were recorded by a thermocouple calibrated in magnetic 
field using measurements on ultrapure silicon. We per- 
formed measurements at temperatures down to 0.6 K and 
in magnetic fields up to 8 T in the 3 He and 4 He refriger- 
ators in Kosice. Supplementary measurements up to 14 
Tesla and down to 2 K were performed in Grenoble. 
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FIG. 1: a) Total specific heat C/T of NbS2 in magnetic field 
parallel (a) and perpendicular (b) to the ab planes. In both 
data sets the zero field measurement is the rightmost curve. 



III. RESULTS AND DISCUSSION 

Figure 1 displays the temperature dependence of the 
specific heat of the sample (plus addenda) in selected 
magnetic fields up to 8 T for H \ \ab and up to 3 T for H\\c. 
The thermodynamic superconducting transition temper- 
ature at zero field was determined from the local entropy 
balance around the phase transition giving T c — 6.05 K. 
The zero-field anomaly at the transition is sharp ( AT C 
~ 0.4 K) indicating the high quality and homogeneity of 
the sample. The position of the specific-heat jumps are 
gradually shifted toward lower temperatures for increas- 
ing magnetic field. Despite a significant broadening at 
high fields, the anomaly remains well resolved at all fields. 
A field of 3 Tesla applied along the c-axis was sufficient 
to completely suppress superconductivity in all the tem- 
perature range. On the other hand, 8 Tesla applied along 
the ab-planes shifts the superconducting anomaly down 
to only about 3-4K underlying the strong anisotropy of 
this system. 

Later we extended the measurements down to 0.6 K 
in a 3 He fridge where the specific heat was measured at 
zero field and at H\ |c = 3 T. In the case of a very small 
crystal like ours, it is difficult to evaluate the exact total 
contribution of the addenda. The electronic part of the 



total specific heat value can be obtained by extrapolation 
oiCtot/T for T approaching zero. This value corresponds 
to ~ 38% of 7„, with j n T being the electronic heat ca- 
pacity of the sample in the normal state. To avoid any 
fitting procedures, the addenda and the phononic con- 
tributions have been eliminated by subtracting the data 
taken at H\\c = 3 T from all the other runs. Thus, 
the electronic specific heat of the sample, C e (T), nor- 
malized to its normal state value, 7«T, can be obtained 
experimentally by: C e (T)/^ n T = A ( c / T ) _|_ i ; where 
A(C/T) = {C{T.H = 0)/T-C{T,H = 3T||c)/T and 7 „ 
= C(H = 3T\\c)/T ]0 _ 6K -C{H = 0)/Tj . 6Jf . The result is 
presented in Fig. 2 by the open circles. The only assump- 
tion in this procedure is the absence of magnetic field 
dependence of the addenda. This has been previously 
verified in numerous experiments using the same ther- 
mocouple wires and also confirmed here independently 
by the entropy conservation required for a second order 
phase transition, proving the thermodynamic consistency 
of the data and its treatment. 

We first compare the electronic specific heat with the 
isotropic single gap (ISG) BCS model. The dashed line 
in Fig. 2 presents the ISG BCS specific heat (weak cou- 
pling of 2A/fcsT c = 3.52). One can see that while the 
height of the jump at T c of the experimental data is quite 
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FIG. 2: Open circles: Electronic specific heat of NbS2 in zero 
magnetic field extended down to 0.6 K. Dashed line: BCS 
single-gap weak-coupling case. Solid line: two-gap model with 
2As/fcflT c = 2.1, 2Az,/fc.B T c = 4.6 and respective relative 
contributions js/jn = 0.4, 7^/7™ = 0.6. The anisotropic gap 
model with anisotropy parameter a — 0.5 and 2Ao/ksT c = 
3.6 follows essentially the same line. Inset: exponential de- 
pendence of the specific heat, the full line represents the best 
fit of the exponential decay, the dashed line is the behavior 
expected for a BCS single-gap weak-coupling limit. 



well reproduced, a significant deviation occurs at lower 
temperatures. The discrepancy between the BCS curve 
and the measured data represents 18 % (7 %) of the total 
signal at I.4K (4K), which is significantly larger than the 
error bars of our measurements. 

The inset of Fig. 2 displays the logarithm of the elec- 
tronic specific heat versus T c /T. As shown, one ob- 
tains an exponential dependence C e oc exp(—b * T c /T) 
for T c /T > 2.5. However, the parameter b is significantly 
lower than the value expected for the BCS weak cou- 
pling limit in the temperature range T c /T — 2.5 — 4.5 
flq . This corresponds to coupling ratio 2A/ksT c ~ 2.3 
that is much smaller than the ISC BCS value of 3.52, in- 
dicating that the quasiparticles are activated over a small 
energy gap. This fact as well as the overall shape of the 
specific heat temperature dependence resembles the case 
of MgB 2 . 

A phenomenological a-model of the specific heat ac- 
counting for independent contributions from two bands 
with two different energy gaps has been successfully ap- 
plied in the case of MgB2 [l9(. The magnitude of the 
small gap 2As//cbT c and of the large gap 2Ai/fcsT c at 
T = are fitting parameters of the model. The third 
parameter is the relative fraction of the density of states 
of the two bands "ys,n/lL.n- The full line in Fig. 2 repre- 



sents a fit to the experimental data yielding the follow- 
ing parameters: 2A s /k B T c = 2.1 ±0.05, 2A L /k B T c = 
4.6±0.2, and 7s,„/7L, n = 0.67± 0.15. The value of the 
small gap is close to the one evaluated from the exponen- 
tial decay shown in the inset of Fig. 2. Importantly, both 
gap values are in striking agreement with those found in 
the STM experiment [15j confirming that the latter are 
characteristic of the bulk. 

As previously shown by Huang et al. [lOj, the tem- 
perature dependence of the specific heat of NbSe2, an- 
other two-gap superconductor, can also be described by 
an anisotropic s-wave model, where the gap anisotropy is 
supposed to be in the form of A = Ao(l + a cos 6^) cor- 
responding to the hexagonal in-plane symmetry. Here, 
Ao is the average gap value and a denotes its anisotropy, 
yielding A max = A (l + a) and A mm = A (l - a). This 
model with parameters a = 0.5 and 2Ao/ksT c — 3.6 fits 
our data as well as the two gap scenario, the dif ference 
between the two models is negligible. We remark that 
the anisotropic gap should leave its footprint in the 
anisotropic vortex core £ as it is proportional to the re- 
lated Fermi velocity vf divided by the gap at zero tem- 
perature A(0). However, in contrast to NbSe2, for which 
STS images revealed a sixfold star shape of the vortex 
cores, the fully isotropic vortices have been imaged in 
NbS 2 [f 5] questioning the applicability of the anisotropic 
single gap model in our case. 

We have measured thoroughly the evolution of the spe- 
cific heat in the mixed state. At T = 0.6 K the electronic 
specific heat term C e /T is very close to the Sommerfeld 
coefficient 7. Its field dependence is displayed in Fig. 3a 
and 3b for both principal field orientations. Our maxi- 
mum field available at this temperature range (8 Tesla) 
was not sufficient to reach the normal state for ff||a&, 
but it was well above the upper critical field value of 
2.4±0.I Tesla for H\\c. Fig. 3a emphasizes the strong 
non-linearity of j(H) when H is applied perpendicular to 
the ab plane. Again, such a non-linearity could be asso- 
ciated with the existence of 2 gaps or a single anisotropic 
one. 

The increase of 7 with magnetic field is mainly due to 
the quasiparticle contribution inside the vortex cores. In 
the case of superconductor with a single isotropic gap, 
7 should increase linearly in small magnetic field and a 
small nonlincarity in "/(H) appears above the field where 
flux lines start overlapping. According to the calcula- 
tions of Nakai et al. [2pJ, much stronger non-linearity 
of j(H) is achieved in case of anisotropic-gap supercon- 
ductors. The full line in the Fig. 3a displays the field 
dependence of the normalized density of states N(B)/Nq 
(proportional to the Sommerfeld coefficient) calculated 
by Nakai et al. for the anisotropic gap with a = 0.5. The 
model qualitatively reproduces the behavior in our data 
(open circles), but fails to describe the fast increase of 
7(H) at low fields. This is evident in the inset of Fig. 3a 
where the derivative dN(B)/dB, i.e. the slope of N(B) 
for a = 0.5 (line) and the slope of measured Sommerfeld 
coefficient d-f(H)/dH (open circles) are compared. This 
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FIG. 3: (a) Open circles - normalized Sommerfeld coefficient 7 
as a function of magnetic field perpendicular to the ab planes 
of NbS2. Line - model accounting for highly anisotropic gap 
with a = 0.5 [2(j. Inset is the derivative of the corresponding 
curves from the main panel: open circles - of the measured 
data, line - of the model, (b) and (c) 7/7^ for both orien- 
tations of the magnetic field in NbS2 and MgEb [2l[, respec- 
tively. 



discrepancy makes the explanation of the specific heat 
data by the anisotropic gap scenario rather inconsistent. 
The observed 7(1?) behavior resembles the two-gap case 
of MgB2 (see Fig. 3c, left curve) [2l|, where the initial 
rapid increase of 7(H) due to dominant role of the tt 
band with the small gap changes at fields where the a 
band with the large gap comes into play. 

Next, we compare the behavior of the Sommerfeld co- 
efficients in magnetic fields applied parallel and perpen- 
dicular to the basal planes of the both materials, NbS2 
(Fig. 3b) and MgB 2 (Fig. 3c). The latter is taken from 




FIG. 4: Anisotropy of NbS2 (full circles) compared to MgB2 
[2lL I27I ] (open circles): (a) field dependence of effective 
anisotropy defined as the ratio of the fields applied in both 
principal orientations that correspond to the same 7 value in 
Fig. 3b and Fig. 3c (b) temperature dependence of anisotropy 

r = 



Ref.[2l|. Following the procedure introduced by Bouquet 
et al. in Ref. [22| an effective anisotropy T e f / can be ob- 
tained from these 7(H) dependences. It is defined as the 
ratio of the magnetic fields in the ab plane and along the 
c axis yielding the same 7 value in Fig. 3b and Fig. 3c. 
T e t f is plotted in Fig 4a for both compounds as a func- 
tion of the field H\\c normalized to its upper critical field 
value H C 2 \\c. As discussed already in the work of Bouquet 
et al. on MgB2 [HJ the choice of the abscissa is arbitrary. 
We chose H\\c, but we could plot T e ff versus H\\ab or 
versus 7(H) as well. Note that this T e ff tends towards 

H ab 

the usual anisotropy of H C 2, T = -jpr- when -2- — > 1 at 

large magnetic fields. In MgB 2 at low fields <-^- < 0.1), 

the 7(H) curves for the two principal directions are prac- 
tically identical which gives F e f f = 1 as shown in Fig. 4a. 
At larger fields, T e ff increases reflecting a reduced con- 
tribution from the isotropic 7r-band, reaching T e f r ~ 5 
which is the anisotropy of the dominant a-band [2l| . In 
NbS 2 , one observes an opposite field dependence of T e f f 
which starts from a highly anisotropic value T e ff ~ 10 at 
low fields and decreases to T e ff ~ 5.5 at our maximum 
field. A field dependent superconducting anisotropy is a 
typical signature of multigap superconductivity where a 
role of bands with different gaps can significantly vary 
with magnetic field [23h25| . In contrast to MgB 2 case, in 
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FIG. 5: Upper critical field H C 2 for magnetic field \\ab (full 
circles result from temperature-sweep measurements, full tri- 
angles from field sweeps) and \\c (open circles result from 
temperature-sweep measurements, open triangles from field 
sweeps) . 



Importantly, the three independent sets of tempera- 
ture and field-sweeps measurements yield consistent re- 
sults with a nice overlap. Both temperature dependen- 
cies show a slight positive curvature for temperatures 
T > T c /2. The upper critical field in the ab plane re- 
veals very high values with dH c ijdT slope of about 3 
Tesla/K, close to the Pauli paramagnetic limit. 

The temperature dependence of the superconducting 
anisotropy T calculated as a ratio H c2 \\ab/H C 2\\c of the 
upper critical fields is displayed in Fig. 4b together with 
the results obtained in MgB2 [27j . As shown, in con- 
trast to MgB2, the resulting anisotropy T is close to 7 
and approximatively constant for T/T c > 0.3. Note that 
this value might be slightly underestimated in case of a 
small misalignement of the crystal for H\\ab. Our re- 
sults are consistent with those obtained by Onabe et al. 
28] from resistive measurements in a field up to 2 Tesla. 
The strong decrease of T in MgB2 close to T c is a direct 
consequence of the existence of the isotropic n- band. 
This is not a general feature of multiband superconduc- 
tivity since T(T) results from a subtle balance between 
the Fermi velocities and the relative weight in the DOS 
of the different bands [13] . These precise calculations are 
still to be carried out in the case of NbS2 . 



NbS2 both bands would be anisotropic, as suggested by 
analogy with NbSc2 in which band structure calculations 
26] show mostly 4 Fermi surface sheets derived from Nb 
d-bands forming warped cylinders along the c axis, cen- 
tered on the r and K points in the Brillouin zone. More- 
over, two sheets derived from the bonding Nb d band are 
significantly more warped than the two derived from the 
antibonding Nb d band. Different warping of Nb sheets 
can naturally lead to a different level of anisotropy in 
each band. Thus, a qualitatively different behavior of 
T e ff(H) compared to MgB 2 can be expected. 

Finally, we inspected the upper critical magnetic 
fields for both principal orientations of magnetic field. 
Fig. 5 summarizes the values of H C 2 derived from the 
temperature-sweep measurements of the specific heat 
shown in Fig. 1, as well as from field-sweep measure- 
ments. Two sets of field-sweeps were performed, one in 
14 Tesla magnet for if||a& in a temperature range down 
to 2 K, and another one taken in the 3 He cryostat down 
to 0.6 K in the 8 Tesla coil. As stated above, we de- 
termined T c at zero field from the local entropy balance 
around the anomaly. However, at finite fields this def- 
inition is not very practical for establishing T C 2(H), or 
H C 2(T). In order to reduce the uncertainty of the H C 2 
value arising from the broadening of the transition partic- 
ularly at lower temperatures (higher fields) , we inspected 
the temperature shift between two neighboring curves in 
Fig. 1. A similar procedure was used to determine H c2 
from field-sweep measurements. The resulting tempera- 
ture dependence of H C 2 is shown in Fig. 5 for both H \ \ab 
and H\\c. 



IV. CONCLUSIONS 

Analysis of the zero-field specific heat data has shown 
that: (1) the zero-field electronic term of the specific heat 
cannot be described by an isotropic single-gap BCS for- 
mula, but it is compatible with the two-gap a model; (2) 
the large (small) gap is 2A L /k B T c m 4.6 (2A s /k B T c « 
2.1). The measurements in the mixed state have sup- 
ported the two-gap scenario revealing: (3) a strongly 
non-linear j(H); (4) a field-dependent superconducting 
anisotropy. Even if some of these features of the spe- 
cific heat could be eventually explained by an extremely 
anisotropic-gap superconducting model, this would not 
be compatible with the observation of two well-resolved 
gap features with sizes of 2As//cbT' c « 2 and 4, respec- 
tively, and also with the absence of in-plane anisotropy 
in the vortex lattice images by the scanning-tunneling 
spectroscopy of Guillamon et al. To conclude, our bulk 
thermodynamic measurements are in full agreement with 
STM spectra, supporting that NbS2 is another case of 
well resolved two-gap superconductor. 
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